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Q I hidden bottom tetraquark states which consist of the axial-axial type and the 
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CIh 1 Introduction 

I) : 

I— li In 2007, a distinct peak (Z(4430)) was observed in the n^ip' invariant mass dis- 

^ I tribution near 4.43 GeV in the decays B KTc^ip' by the Belle collaboration 

> ; [T]. The fitted Breit-Wigner mass and width are Mz = 4433 ± 4 ± 2MeV and 

^ ■ Tz = 4:5^.11^.13 MeV. The statistical significance of the observed peak is 6.5 a. Using 

Psj . the same data sample, the Belle collaboration also performed a full Dalitz plot anal- 

^ J ! ysis with a fitted model that takes into account all the known Kit resonances below 

^ I 1780 MeV [2]. The significance of the fitted resonance is of 6.4 cr and agrees with 

' the previous observation pQ, the updated parameters are Mz = (44431^2-13) MeV 

O ■ and Tz = {109lill~^lQ) MeV. However, the BaBar collaboration do not confirm this 

. resonance [3], i.e. they observe no significant evidence for a Z(4430) signal for any 

^ ! of the processes investigated, neither in the total J/ ipTc or ip'-K mass distribution nor 

^ I in the corresponding distributions for the regions of Kti mass for which observa- 

■ - - ' tion of the Z(4430) signal is reported. If the Z(4430) exists indeed, it can't be a 

pure cc state due to the positive charge, and may be an excellent tetraquark state 
{ccud) candidate [U [5] . We can distinguish the multiquark states from the hybrids 
or charmonia with the criterion of non-zero charge. 

In 2008, the Belle collaboration reported the first observation of two resonance- 
like structures (the Z(4050) and Z(4250)) in the vr+Xci invariant mass distribution 
near 4.1 GeV in the exclusive decays K~n^Xci [Sj. Their quark contents must 

be some special combinations of the ccud, just like the Z(4430), they can't be the 
conventional mesons. The Z(4050) and Z(4250) lie about (0.5 — 0.6) GeV above 
the vr+Xci threshold, the decay Z t^^Xcv can take place with the "fall-apart" 
mechanism and it is OZI super-allowed, which can take into account the large total 
width naturally. The spins of the Z(4050) and Z(4250) are not determined yet, 
they can be scalar or vector mesons. If they are scalar mesons, the decays Z 
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7T~^Xci occur through the relative P-wave with the phenomenological lagrangian C = 
gx'^iT^daZ — ZdaTi). On the other hand, if they are vector mesons, the decays occur 
through the relative S-wave with the phenomenological lagrangian C = gx°'Za7r. 
There have been several interpretations, such as the tetraquark states [3 El [9] and 
the molecular states [TOl [TH [T2| [13] . 

In Refs.[7l[8], we assume that the hidden charm mesons Z(4050) and Z(4250) 
are vector (and scalar) tetraquark states, and study their masses with the QCD 
sum rules. The numerical results indicate that the mass of the vector hidden charm 
tetraquark state is about Mz = (5.12 ± 0.15) GeV or Mz = (5.16 ± 0.16) GeV, 
and the mass of the scalar hidden charm tetraquark state is about Mz = (4.36 ± 
0.18) GeV. In Refs.fHl |15], we study the mass spectrum of the scalar and vector 
hidden charm and hidden bottom tetraquark states using the QCD sum rules, and 
observe that the scalar hidden charm tetraquark states may have smaller masses 
than the corresponding vector states. From our previous works, we can draw the 
conclusion that the hidden charm meson Z(4250) may be a scalar tetraquark state 
[3 El [HI [15], although other possibilities, such as a hadro-charmonium resonance 
and a D^D^ + D^D^ molecular state are not excluded. We intend to study the mass 
spectrum of the scalar hidden charm and hidden bottom tetraquark states which 
consist of diquark pairs differ from our previous works. 

The mass is a fundamental parameter in describing a hadron, whether or not 
there exist those hidden charm and hidden bottom tetraquark configurations is of 
great importance itself, because it provides a new opportunity for a deeper under- 
standing of the low energy QCD. 

In Refs. p!B| [T7] , Ebert et al take the diquarks as bound states of the light and 
heavy quarks in the color antitriplet channel, and calculate their mass spectrum 
using a Schrodinger type equation, then take the masses of the diquarks as the ba- 
sic input parameters, and study the mass spectrum of the heavy tetraquark states 
as bound states of the diquark-antidiquark system. In Refs.[TSl [T^l I2D], Maiani et 
al take the diquarks as the basic constituents, examine the rich spectrum of the 
diquark-antidiquark states with the constituent diquark masses and the spin-spin 
interactions, and try to accommodate some of the newly observed charmonium- 
like resonances not fitting a pure cc assignment. In Ref. [21j. Zouzou et al solve 
the four-body (QQqq) problem by three different variational methods with a non- 
relativistic potential considering explicitly virtual meson-meson components in the 
wave-functions, search for possible bound states below the threshold for the sponta- 
neous dissociation into two mesons, and observe that the exotic bound states QQqq 
maybe exist for unequal quark masses (the ratio niQ/mq is large enough). The stud- 
ies using a potential derived from the MIT bag model in the Born-Oppenheimer 
approximation support this observation [221 E3]- In Ref.pi], Manohar and Wise 
study systems of two heavy-light mesons interacting through an one-pion exchange 
potential determined by the heavy meson chiral perturbation theory and observe 
the long range potential maybe sufficiently attractive to produce a weakly bound 
two- meson state in the case Q = b. In Ref.|2S]) th^ L = tetraquark states QQQQ 
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{Q denotes both Q and q) are analyzed in a chromo- magnetic model where only a 
constant hyperfine potential is retained. 

In this article, we re-study the mass spectrum of the scalar hidden charm and 
hidden bottom tetraquark states using the QCD sum rules [26l|27j. In the QCD sum 
rules, the operator product expansion is used to expand the time-ordered currents 
into a series of quark and gluon condensates which parameterize the long distance 
properties of the QCD vacuum. Based on the quark-hadron duality, we can obtain 
copious information about the hadronic parameters at the phenomenological side 

The hidden charm and hidden bottom tetraquark states (Z) have the symbolic 
quark structures: 

Z+ = QQud; Z° = -^QQiuu - dd); Z' = QQdu; 

= QQus; Z~ = QQsu; Z^ = QQds; Z ^ = QQsd; 
Z^ = -^QQiuu + dd)] Z^ = QQss, (1) 

where the Q denotes the heavy quarks c and b. 

We take the diquarks as the basic constituents to study the tetraquark states fol- 
lowing Jaffe and Wilczek ES] ■ The heavy tetraquark system could be described 
by a double-well potential with two light quarks q'q lying in the two wells respec- 
tively. In the heavy quark limit, the c (and b) quark can be taken as a static well 
potential, which binds the light quark q to form a diquark in the color antitriplet 
channel. The attractive interactions of one-gluon exchange favor formation of the 
diquarks in color antitriplet 3c, flavor antitriplet 3/ and spin singlet Ig [301 El]. The 
diquarks have five Dirac tensor structures, scalar C75, pseudoscalar C, vector C'j^'j^, 
axial vector 6*7^ and tensor Ccr^i,. The structures 6*7^ and Ca^u are symmetric, 
the structures C75, C and C'^f^'y^ are antisymmetric. In Refs. [H HI] , we assume the 
scalar hidden charm and hidden bottom mesons Z consist of the C75 — C75 type 
diquark structures rather than the C — C type diquark structures, and observe that 
the C75 — C75 type tetraquark states have much smaller masses than the corre- 
sponding C — C type tetraquark states; our numerical results of the C — C type 
tetraquark states will be presented elsewhere. 

In this article, we assume the scalar hidden charm and hidden bottom tetraquak 
states which consist of the 6*7^ — C'j'^ type and the 6*7^75 — Cj'^j^ type diqaurk 
pairs and study the mass spectrum. Naively, we expect the 6*7^ — 6*7'^ type and the 
C7At75 — C*7^75 type tetraquark states have larger masses than the corresponding 
C75 — C75 type tetraquark states. 

The article is arranged as follows: we derive the QCD sum rules for the scalar 
hidden charm and hidden bottom tetraquark states Z in section 2; in section 3, 
numerical results and discussions; section 4 is reserved for conclusion. 
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2 QCD sum rules for the scalar tetraquark states 
Z 

In the following, we write down the two-point correlation functions n(p) in the QCD 
sum rules, 

n(p) = i J d^xe'P^O\T{j/ri{x)J/ri\0)}\0) , (2) 

where the J{x) and ri{x) denotes the interpolating currents J^+(x), Jzo{x), 
Vz+{x), Vzo{x), • • •, 

Jz+{x) = e'^h''^^uJ{x)C^^Qk{x)QUx)YCd:^{x), 

Jzo{x) = ^ [uj {x)C'j^Qkix)Qmix)YCul{x) -{u^ d)] , 

JzA^) = e'"'e'"'^uJ{x)Cj,Qk{x)QUxh''Csl{x), 
Jzo{x) = e''''e'^^dJ{x)C^M^)Qm{xh^Csl{x), 

Jz^x) = ^ [uj {x)C-f^Qk{x)Qmix)YCul{x) + (m d)] , 

Jz,{x) = e'^'e'"'^sJ{x)Cj,Qk{x)QUxh''Csl{x), 
Vz+{x) = e'"'e''^^uJ{x)C^,-f,Qk{x)QUxh5YCd^{x), 

^ijk^imn 

VZ»{x) = ^ [uJ{x)C-ff,-f5Qk{x)Qrn{x)-f5l''Cul{x) - {u ^ d)] , 

r7^+(x) = e'"'e'"'^uJ{x)Cj,j,Qk{x)QUxHl''Csl{x), 
r^zo{x) = e'^''e'"'^dJ{x)C^,j,Qk{x)QUxHYCsl{x), 

^ijk^imn 

Tlz^{x) = ^ [uJ{x)C-ff,-f5Qk{x)Qm{x)-f5l''Cul{x) + {u d)] , 

Vz,{x) = e''h'"'''sJ{x)C^,^,Qk{x)QUxHYCsl{x), (3) 

where the i, j, k, ■ ■ ■ are color indexes. In the isospin limit, the interpolating 
currents result in six distinct expressions for the correlation functions 11 (p) , which 
are characterized by the number of the s quark they contain. In Refs. [T3] . 
we observe that the ground state masses of the scalar and vector tetraquarks are 
characterized by the number of the s quarks they contain, Mq < Mg < Mgs] the 
energy gap between Mq and Mss is about (0.05 — 0.15) GeV. In this article, we study 
the interpolating currents which contains zero and two s quarks for simplicity. 

We can insert a complete set of intermediate hadronic states with the same quan- 
tum numbers as the current operators J{x) and ri{x) into the correlation functions 
n(p) to obtain the hadronic representation [261127] . After isolating the ground state 
contribution from the pole terms of the Z, we get the following result, 
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where the pole residue (or couphng) \z is defined by 

\z = (0|J/77(0)|Z(p)). (5) 

The contributions from the two-particle and many-particle reducible states are 
supposed to be small enough to be neglected safely, for example, the scattering state 
XciTi"'' in the ccdu channel. 



' ^^^-^ J (27r)4 [q^ - mlJ [{p - - ml\ 



9lJ.i> ~r 2 



+ •••,(6) 



where 

(0|Jz+(0)|Xci7r+) = A^^,,+p^e^ (7) 

the is the polarization vector of the axial- vector meson Xci- We can estimate the 
coupling A^^j7r+ with the soft tt theorem, 

(0|Jz+(0)|xci7r+) = -f (0|[g5,Jz+(0)]|Xci), 

/tt 

= f(0|Jp(0)|Xci) = -^, (8) 

/tt /tt 

where 

Jp{x) = e'^h''^^[u]{x)C-i^-f^Ck{x)crn{xh''Cul{x) 

+4(x)C7^Cfe(x)c^(x)757^C(i^(^)] . (9) 

As the main Fock states of the charmonia are the cc components, the coupling Ap 
between the pseudoscalar tetraquark current Jp{x) and the axial- vector meson Xci 
should be very small. 

We can perform Fierz re-ordering in both the Dirac spin space and the color 
space to express the tetraquark current Jz+ {x) in the following form, 

3 - 3 - 3 - 

Jz+{x) = -d{x)u{x)c{x)c{x) -d{x)i'y^u{x)c{x)i'y^c{x) —d{x)-fau{x)c{x)-i°'c{x) 
8 8 Id 

3 - 1 - A* A* 

-YgO?(2^)7a75^iWc(a;)7"75c(a;) - -d{x)—u{x)c{x)—c{x) 

1 - A* A* 1 - A* A* 

--d(x)i75— M(a;)c(x)i75— c(x) - -d{x)-^a—u{x)c{x)-^'^—c{x) 

1 - A* A* 

+-d(a:)7a75y?i(x)c(x)7"75— c(x) , (10) 

where the A* are the matrix elements of the SU{2) group in adjoint representation, 
the i = 1,2,3 are the color indexes. The scalar tetraquark current which consists 
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of a axial-vector diquark pair is a special composition of the S — S,P — P,V — V, 
A — A, — S\ — P\ — V^ and A^ — A' color-singlet and color-triplet meson- 
meson type currents, the S, P, V and A denote the scalar, pseudoscalar, vector and 
axial-vector respectively. The color-singlet meson-meson type currents, for exam- 
ple, (J(x)7q,75m(x)c(x)7"75c(x), J(x)i75'u(x)c(x)275c(x), have very small two-particle 
reducible contributions [32] . 

After performing the standard procedure of the QCD sum rules, we obtain the 
following four sum rules for the interpolating currents contain two s quarks: 



the explicit expressions of the spectral densities p±{s) are presented in the appendix, 
the -|- and — denote the Cj^ — C'j'^ type and the 6*7^75 — C'j'^j^ type interpolating 
currents respectively; the Sq is the continuum threshold parameter and the is 
the Borel parameter. We can obtain four sum rules in the ccqq and bbqq channels 
with a simple replacement rrig, (ss) — > {qq) and {sggCrGs) {qgsCrGq). 

We carry out the operator product expansion to the vacuum condensates adding 
up to dimension-10. In calculation, we take assumption of vacuum saturation for 
high dimension vacuum condensates, they are always factorized to lower condensates 
with vacuum saturation in the QCD sum rules, factorization works well in large Nc 
limit. In this article, we take into account the contributions from the quark conden- 
sates, mixed condensates, and neglect the contributions from the gluon condensate. 
The contributions from the gluon condensates are suppressed by large denominators 
and would not play any significant roles for the light tetraquark states [331 El] , the 
heavy tetraquark state [8j and the heavy molecular state [35] . There are many terms 
involving the gluon condensate for the heavy tetraquark states and heavy molecu- 
lar states in the operator product expansion (one can consult Refs.[Sl ESI ESj for 
example), we neglect the gluon condensates for simplicity. 

Differentiate the Eq.(ll) with respect to then ehminate the pole residues 
Xz, we can obtain the sum rules for the masses of the tetraquark quark states Z, 



3 Numerical results and discussions 

The input parameters are taken to be the standard values (qq) = — (0.24±0.01 GeV)'^, 
{ss) = (0.8 ± 0.2) {qq), {qg^crGq) = ml{qq), {sg.aGs) = ml{ss), ml = (0.8 ± 
0.2) GeV^ m, = (0.14 ± 0.01) GeV, = rrid ^ 0, = (1.35 ± 0.10) GeV and 
nib = (4.8 ± 0.1) GeV at the energy scale /i = 1 GeV [26l [271 [37] . 

In the conventional QCD sum rules |26l [27|, there are two criteria (pole dom- 
inance and convergence of the operator product expansion) for choosing the Borel 




(11) 



j2 dsp±{s)e'^ 



(12) 
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parameter and threshold parameter sq. We impose the two criteria on the heavy 
tetraquark states to choose the Borel parameter and threshold parameter sq. 

In Refs.[5l[n]. we assume that the resonance-like structures ^(4050) and Z(4250) 
are scalar tetraquark states which consist of the scalar diquark pairs, and take the 
threshold parameter tentatively as sq = (4.248 + 0.5)^ GeV^ ^ 23 GeV^ to take 
into account all possible contributions from the ground states, where the energy gap 
between the ground states and the first radial excited states is chosen to be 0.5 GeV. 
Then we take into account the SU{3) symmetry of the light flavor quarks and 
the mass difference between the heavy quarks, choose other threshold parameters 
tentatively, and use those values as a guide to determine the threshold parameters 
So with the QCD sum rules. 

In this article, we study the scalar hidden charm and hidden bottom tetraquark 
states which consist of the axial-axial type and the vector-vector type diquark 
pairs, and search for other possible tetraquark structures of the resonance- like states 
Z(4050) and Z(4250). Naively, we expect the C-ff,-f5-CY75 type and the C-ff,-CY 
type tetraquark states have larger masses than the corresponding C75 — C75 type 
tetraquark states, and use the threshold parameters in Ref. [Hj as a guide to deter- 
mine the threshold parameters Sq with the QCD sum rules. 

The contributions from the high dimension vacuum condensates in the operator 
product expansion are shown in Figs. 1-2, where (and thereafter) we use the (qq) 
to denote the quark condensates {qq), {ss) and the {qggCrGq) to denote the mixed 
condensates {qggaGq), {sggCrGs). From the figures, we can see that the contributions 
from the high dimension condensates change quickly with variation of the Borel 
parameter at the values < 2.6 GeV^ (2.8 GeV^) and < 7.2 GeV^ (7.6 GeV^) 
in the hidden charm and hidden bottom channels respectively for the 6*7^ — 6*7'^ 
(6*7^75 — C7'^75) type interpolating currents, such an unstable behavior cannot lead 
to stable sum rules, our numerical results confirm this conjecture, see Fig. 4. 

At the values > 2.6 GeV^ (2.8 GeV^) and Sq > 23 GeV^ (27GeV^), the con- 
tributions from the (gg)^ + {qq){qgs<^Gq) term are less than 14% (23.5%) in the 
ccqq channel, the corresponding contributions are less than 4% (13%) in the cess 
channels; the contributions from the vacuum condensate of the highest dimension 
{qgsO'Gq)'^ are less than 2.5% (2.5%) and 1.5% (3%) in the ccqq and cess channels re- 
spectively; we expect the operator product expansion is convergent for the C-y^—C'j^ 
(C*7At75 ~ C'-f'^'-f^) type interpolating currents in the hidden charm channels,. 

At the values > 7.2 GeV^ (7.6 GeV^) and sq > 136 GeV^ (146 GeV^), the 
contributions from the (qq)"^ + {qq){qgsO'Gq) term are less than 10.5% (19%) in the 
bbqq channel, the corresponding contributions are less than 3.5% (8.5%) in the bbss 
channels; the contributions from the vacuum condensate of the highest dimension 
{qgscGqy are less than 5% (6%) and 3% (6%) in the bbqq and bbss channels respec- 
tively; we also expect the operator product expansion is convergent for the C'jf^—C'y'^ 
(C'TajTs — ^7^75) type interpolating currents in the hidden bottom channels. 

In this article, we take the uniform Borel parameter M^^^, i.e. M^j„ > 2.6 GeV^ 
(2.8 GeV^) and M^^^ > 7.2 GeV^ (7.6 GeV^) in the hidden charm and hidden bottom 
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channels respectively for the C'-f^—Cj'^ (CT/xTs^Ct'^Ts) type interpolating currents. 

In Fig. 3, we show the contributions from the pole terms with variation of the 
Borel parameters and the threshold parameters. The pole contributions are larger 
than (or equal) 50% (52%) at the value < 3.2 GeV^ and sq > 23 GeV^ (27 GeV^), 
24 GeV^ (28 GeV^) in the ccqq, cess channels respectively, and larger than (or equal) 
51% (52%) at the value M\< 8.2 GeV^ and Sq > 136 GeV^ (146 GeV^), 138 GeV^ 
(148 GeV^) in the bbqq and bbss channels respectively for the the C'j^—C'j^ (C7/^75~ 
^7^75) type interpolating currents. Again we take the uniform Borel parameter 
M^„^, i.e. M^„^ < 3.2 GeV^ and M^^^ < 8.2 GeV^ in the hidden charm and hidden 
bottom channels respectively. 

For the 6*7^ — 6*7^ type interpolating currents, the threshold parameters are 
taken as Sq = (24 ± 1) GeV^ (25 ± 1) GeV^ (138 ± 2) GeV^ and (140 ± 2) GeV^ in 
the ccqq, cess, bbqq, and bbss channels respectively; the Borel parameters are taken 
as = (2.6 - 3.2) GeV^ and (7.2 - 8.2) GeV^ in the hidden charm and hidden 
bottom channels respectively. 

For the 6*7^75 — 6*7^^75 type interpolating currents, the threshold parameters are 
taken as Sq = (28 ± 1) GeV^ (29 ± 1) GeV^ (148 ± 2) GeV^ and (150 ± 2) GeV^ for 
the ccqq, cess, bbqq, and bbss channels, respectively; the Borel parameters are taken 
as = (2.8 - 3.2) GeV^ and (7.6 - 8.2) GeV^ in the hidden charm and hidden 
bottom channels respectively. 

In those regions, the pole contributions are about (47 — 75)%, (51 — 78)%, (51 — 
70)% and (53 — 72)% in the ccqq, cess, bbqq and bbss channels respectively for 
the 6*7^ — C'-y^ type interpolating currents; while the pole contributions are about 
(52 - 75)%, (52 - 74)%, (52 - 68)% and (52 - 67)% in the ccqq, cess, bbqq and bbss 
channels respectively for 6*7^75 — 6*7^75 type interpolating currents; the two criteria 
of the QCD sum rules are fully satisfied [26i I27j . 

If we take uniform pole contributions, the interpolating current with more s 
quarks requires slightly larger threshold parameter due to the SU{3) breaking effects, 
see Fig. 3. The threshold parameters in the ccqq and bbqq channels are slightly 
smaller than the corresponding ones in the cess and bbss channels respectively. 
Naively, we expect the tetraquark state with more s quarks will have larger mass, 
our numerical calculations confirm this conjecture, see Fig. 4. In that figure we plot 
the tetraquark state masses Mz with variation of the Borel parameters and the 
threshold parameters. 

The Borel windows M^^^, — M^j„ change with variations of the threshold parame- 
ters So, see Fig. 3. In this article, the Borel windows are taken as 0.6 GeV^ (0.4 GeV^) 
and 1.0 GeV^ (0.6 GeV^) in the hidden charm and hidden bottom channels respec- 
tively for the 6*7^ — 6*7^ (^7/^75 — C7'^75) type interpolating currents; they are small 
enough. Furthermore, we take uniform Borel windows and smear the dependence on 
the threshold parameters sq in each channel. If we take larger threshold parameters, 
the Borel windows are larger and the resulting masses are larger, see Fig. 4. In this 
article, we intend to calculate the possibly lowest masses which are supposed to be 
the ground state masses by imposing the two criteria of the QCD sum rules. 
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tetraquark states 


C7m75 - C'7''75 


c^, - cr 


C75 - C75 


Refs.[T6lfT7] 


cess 


4.82 ±0.14 


4.45 ±0.16 


4.44 ±0.16 


4.110 


ccqq 


4.56 ±0.14 


4.36 ±0.18 


4.37 ±0.18 


3.852 


bbss 


11.70 ±0.18 


11.23 ±0.16 


11.31 ±0.16 


11.133 


bbqq 


11.38 ±0.13 


11.14±0.19 


11.27 ±0.20 


10.942 



Table 1: The masses (in unit of GeV) of the scalar tetraquark states, the values for 
the C75 — C75 type scalar tetraquark states are taken from Ref . [H] . 



tetraquark states 


C7/.75 - C'7''75 


C7m - CY 


cess 


7.92 ± 1.95 


7.05 ± 1.45 


ccqq 


6.32 ±2.30 


5.85 ± 1.30 


bbss 


4.46 ± 1.04 


3.68 ±0.80 


bbqq 


3.35 ± 1.00 


3.06 ±0.66 



Table 2: The pole residues (in unit of 10 ^ GeV^ and 10 ^ GeV^ for the hidden 
charm and hidden bottom channels respectively) of the scalar tetraquark states. 

Taking into account all uncertainties of the input parameters, finally we obtain 
the values of the masses and pole resides of the scalar tetraquark states Z, which 
are shown in Figs. 5-6 and Tables 1-2. In this article, we calculate the uncertainties 
6 with the formula 



where the / denote the hadron mass Mz and the pole residue Xz, the Xi denote 
the input QCD parameters rric, m^, (qq), (ss), As the partial derivatives ^ 

are difficult to carry out analytically, we take the approximation (^§§-^ i^i ~ 

[/(xj ± Axj) — /(xj)]^ in the numerical calculations. 

From Tables 1-2, we can see that the uncertainties of the masses Mz are rather 
small (about (3 — 4)% in the hidden charm channels and (1 — 2)% in the hidden 
bottom channels), while the uncertainties of the pole residues Xz are rather large 
(about (20 — 40)%). The uncertainties of the input parameters {{qq), {ss), {sgsCrGs), 
{qgsCrGq), m^, nic and nib) vary in the range (2 — 25)%, the uncertainties of the pole 
residues Xz are reasonable. We obtain the squared masses M| through a fraction, 
see Eq.(12), the uncertainties in the numerator and denominator which origin from 
a given input parameter (for example, (ss), {sggCrGs)) cancel out with each other, 
and result in small net uncertainty. 

The SU{3) breaking effects for the masses of the hidden charm and hidden 
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bottom tetraquark states are buried in the uncertainties. The — C'y^ type and 
the C75 — C75 type interpolating currents result in almost the same ground state 
masses, while the ground masses of the 6*7^75 — 6*7^75 type tetraquark states are 
larger than the corresponding ones of the C7^ — C'j'^ type tetraquark states about 
(0.2 — 0.5) GeV. Naively, we expect the axial and vector diquarks have larger masses 
than the corresponding scalar diqaurks, and the Cj^ — C'j^ type and the Cj^'j^ — 
C7'^75 type scalar tetraquark states have larger masses than the corresponding C75 — 
C75 type scalar tetraquark states, because the attractive interactions of one-gluon 
exchange favor formation of the diquarks in color antitriplet 3c, flavor antitriplet 3/ 
and spin singlet 1^ [30l |3T] . 

The meson Z(4250) may be a scalar tetraquark state (ccud), irrespective of the 
C'-y^ — 6*7^ type and the C75 — C75 type diquark structures [5], [H], the decay 
Z(4250) TT^Xci can take place with the OZI super-allowed "fall-apart" mecha- 
nism, which can take into account the large total width naturally. Other possibilities, 
such as a hadro-charmonium resonance and a DfD^ + D^D^ molecular state are 
not excluded; more experimental data are still needed to identify it. It is difficult 
to identify the Z(4050) as the scalar tetraquark state as the lower bound of the 
6*7^ — C7^ type and the C75 — C75 type scalar tetraquark states are larger than 
the Z(4050) about 130 MeV. 

In this article, we calculate the mass spectrum of the scalar hidden charm 
and hidden bottom tetraquark states consist of the 6*7^75 — C'j^'y^ type and the 
6*7^ — C'j^ type diquark pairs by imposing the two criteria of the QCD sum rules. 
In fact, we usually consult the experimental data in choosing the Borel parame- 
ter and the threshold parameter sq [5B]- There lack experimental data for the 
phenomenological hadronic spectral densities of the tetraquark states, the present 
predictions can't be confronted with the experimental data. The nonet scalar 
mesons below 1 GeV (the /o(980) and ao(980) especially) are good candidates for 
the tetraquark states [221 EH SO]- However, they can't satisfy the two criteria of 
the QCD sum rules, and result in a reasonable Borel window, although it is not 
an indication non-existence of the light tetraquark states (For detailed discussions 
about this subject, one can consult Refs.[8|i^). The QCD sum rules is just a QCD 
model. 

For the conventional mesons and baryons, the Borel window M^^^ — M^^„ is 
rather large and reliable QCD sum rules can be obtained. However, for the mul- 
tiquark states i.e. tetraquark states, pentaquark states, hexaquark states, etc, the 
spectral densities p ~ s" with n is larger than the ones for the conventional hadrons, 
integral s"' exp{ — ^}ds converges more slowly, which results in some sensitivi- 
ties to the threshold parameters inevitably. We select the threshold parameters and 
Borel parameters by imposing the two criteria of the QCD sum rules, and intend 
to select the possibly lowest threshold parameter which corresponds to the ground 
state. 

In Table 1, we also present the results for the 6*7^ — C'j'^ type scalar tetraquark 
states from a relativistic quark model based on a quasipotential approach in QCD 
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p!6l [T7] , the central values of our predictions are larger than the corresponding ones 
from the quasipotential model about (0.1 — 0.5) GeV. The predications based on 
constituent diquark model {M^cqq = 3723 MeV [IE] and M^sss = 3834 MeV P] for 
the tetraquark states ccqq and cess respectively) are about 0.6 GeV smaller than 
the corresponding ones in the present work. 

The predictions of Refs. [18], [191 [20] depend heavily on the assumption that the 
light scalar mesons ao(980) and /o(980) are tetraquark states, the basic param- 
eters (constituent diquark masses) are estimated thereafter. In the conventional 
quark models, the constituent quark masses are taken as the basic input parame- 
ters, and fitted to reproduce the mass spectra of the well known mesons and baryons. 
However, the present experimental knowledge about the phenomenological hadronic 
spectral densities of the tetraquark states is rather vague, whether or not there ex- 
ist tetraquark states is not confirmed with confidence. The predicted constituent 
diquark masses cannot be confronted with the experimental data. 

The LHCb is a dedicated b and c-physics precision experiment at the LHC (large 
hadron collider). The LHC will be the world's most copious source of the b hadrons, 
and a complete spectrum of the b hadrons will be available through gluon fusion. 
In proton-proton collisions at a/s = 14 TeV, the bb cross section is expected to be ~ 
500/i6 producing 10^^ bb pairs in a standard year of running at the LHCb operational 
luminosity of 2 x 10^^cm~^sec~^ [12]. The scalar tetraquark states (irrespective of 
the 6*7^ — C'-y^ type, the 6*7^75 — C'j^'y^ type and the C75 — C75 type diquark 
structures) may be observed at the LHCb, if they exist indeed. We can search for 
the scalar hidden charm tetraquark states in the DD, D*D*, DsDs, D*D*, J/ipP, 
J/ipip-, J/ipoJ, rjcir, rjcf], ■■■ invariant mass distributions and search for the scalar 
hidden bottom tetraquark states in the BB, B*B*, B^Bs, B*B*, Tp, T0, Tu, rjbTi, 
rjfjT], ■ ■ ■ invariant mass distributions. 

4 Conclusion 

In this article, we study the mass spectrum of the scalar hidden charm and hidden 
bottom tetraquark states which consist of the axial-axial type and the vector- vector 
type diquark pairs with the QCD sum rules, and observe that the scalar-scalar 
type and the axial-axial type tetraquark states have almost the same ground state 
masses while the vector-vector type tetraquark states have slightly larger ground 
state masses. Furthermore, we compare the present predictions with the corre- 
sponding ones from a relativistic quark model based on a quasipotential approach 
in QCD, and discuss the values from the constituent diquark model based on the 
constituent diquark masses and the spin-spin interactions. We can search for the 
scalar hidden charm and bottom tetraquark states at the LHCb. 

We can perform Fierz re-ordering in both the Dirac spin space and the color 
space to express the tetraquark currents J(x) and ri{x) into a series of S — S, P — P, 
V -V, A-A, S'- S\ P' - P\ V' - V' and A' - A' color-singlet and color-triplet 
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Figure 1: The contributions from the high dimension vacuum condensates with 
variation of the Borel parameter in the operator product expansion for the 
6*7^ — 6*7^ type interpolating currents. The (I) and (II) denote the contributions 
from the {qgsCrGq)'^ and (gg)^ +{<l(l){<19sO'Gq) terms respectively. The A, B, C and 
D correspond to the ccgg, cess, bbqq and bbss channels respectively. The notations 
a, /3, 7, A, p and r correspond to the threshold parameters sq = 21 GeV^, 22 GeV^, 
23 GeV^, 24 GeV^, 25 GeV^ and 26 GeV^ respectively in the hidden charm channels; 
while in the hidden bottom channels they correspond to the threshold parameters 
So = 132 GeV^, 134 GeV^, 136 GeV^, 138 GeV^, 140 GeV^ and 142 GeV^ respectively. 
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Figure 2: The contributions from the high dimension vacuum condensates with 
variation of the Borel parameter in the operator product expansion for the 
C'y^'j^ — C'y^'j^ type interpolating currents. The (I) and (II) denote the contributions 
from the {ggscGq)"^ and (gg)^ +{<l(l){Q9sO'Gq) terms respectively. The A, B, C and 
D correspond to the ccgg, cess, bbqq and bbss channels respectively. The notations 
a, /3, 7, A, p and r correspond to the threshold parameters sq = 25 GeV^, 26 GeV^, 
27GeV^, 28 GeV^, 29 GeV^ and 30 GeV^ respectively in the hidden charm channels; 
while in the hidden bottom channels they correspond to the threshold parameters 
So = 142 GeV^, 144 GeV^, 146 GeV^, 148 GeV^, 150 GeV^ and 152 GeV^ respectively. 
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Figure 3: The contributions from the pole terms with variation of the Borel pa- 
rameter M^. The A, B, C, and D denote the ccqq, cess, bbqq, and bbss channels 
respectively. The (I) and (II) denote the 6*7^ — 6*7^ type and the 6*7^75 — C'y^'j^ 
type interpolating currents respectively. For the C'j^—C'y^ (C*7^t75— (^7^75) type in- 
terpolating currents, in the hidden charm channels, the notations a, f3, 7, A, p and r 
correspond to the threshold parameters sq = 21 GeV^ (25 GeV^), 22 GeV^ (26 GeV^), 
23 GeV^ (27 GeV^) , 24 GeV^ (28 GeV^) , 25 GeV^ (29 GeV^) and 26 GeV^ (30 GeV^) , 
respectively; while in the hidden bottom channels they correspond to the threshold 
parameters Sq = 132 GeV^ (142 GeV^), 134 GeV^ (144 GeV^), 136 GeV^ (146 GeV^), 
138 GeV^ (148 GeV^), 140 GeV^ (150 GeV^) and 142 GeV^ (152 GeV^), respectively. 
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Figure 4: The masses of the scalar tetraquark states with variation of the 
Borel parameter and threshold parameter sq. The A, B, C, and D de- 
note the ccqq, cess, bbqq, and bbss channels respectively. The (I) and (II) 
denote the 6*7^ — 6*7^ type and the 6*7^75 — C'-y'^'-y^ type interpolating cur- 
rents respectively. For the 6*7^ — 6*7'^ ((^7^75 — 6*7^75) type interpolating cur- 
rents, in the hidden charm channels, the notations a, (3, 7, A, p and r corre- 
spond to the threshold parameters sq = 21 GeV^ (25 GeV^), 22 GeV^ (26 GeV^), 
23 GeV^ (27 GeV^), 24 GeV^ (28 GeV^), 25 GeV^ (29 GeV^) and 26 GeV^ (30 GeV^), 
respectively; while in the hidden bottom channels they correspond to the threshold 
parameters sq = 132 GeV^ (142 GeV^), 134 GeV^ (144 GeV^), 136 GeV^ (146 GeV^), 
138 GeV^ (148 GeV^), 140 GeV^ (150 GeV^) and 142 GeV^ (152 GeV^), respectively. 



15 




5.5 r 

5.4 ■ 



51 

^ ^ 5.0 - 

> 4.9 

irl,4.7 ■ 

N4,6 : 
^ 4.5 ■ 



B 



Central value; 

Upper bound; 

(I) Lower bound. 



4.3 ■ 
4.2 ■ 



2.0 2.2 2.4 2.6 2.8 3.0 3.2 3.4 3.6 3.8 4.0 



M [GeV ] 




2.2 2.4 2.6 2.8 3.0 3.2 3.4 3.6 3.8 4.0 4.2 



M'^ [GeVT 



12.2 

12.1 
12.0 
11.9 
11.8 
^11.7 
> 11.6 

"ins 
iiii.4 

Nil. 3 
^ 11.2 
11.1 
11.0 
10.9 
10.8 
10.7 









c 


Central value; 






Upper bound; 




(1) 


Lower bound. 



























6.0 6.4 6.8 7.2 7.6 8.0 8.4 

[GeV^] 



9.2 9.6 10.0 



12.0 

11.9 - 

11.8 - 

117 - 



- Central value; C 

- Upper bound; 
Lower bound. (") 




72 7.6 8.0 8.4 8.8 9.2 9.6 10.0 10.4 

[GeV^] 



12.2 
12.1 
12.0 
11.9 
11.8 
117 
^11.6 

CD 

iiii.4 

N11.3 
2 11.2 
11.1 
11.0 
10.9 
10.8 
10.7 







D 


Central value; 




Upper bound; 




(') 


Lower bound. 























7.2 7.6 8.0 8.4 8.8 9.2 

[GeV^] 



12.9 
12.8 
12.7 
12.6 
12.5 
,—,12.4 

(J 12.2 
' — '12.1 
^•^2.0 
11.9 

11.8 
117 
11.6 
11.5 
11.4 



D 
(II) 



-Central value; 
- Upper bound; 
Lower bound. 



7.2 7.6 8.0 8.4 8.8 9.2 9.6 10.0 10.4 

[GeV^] 



Figure 5: The masses of the scalar tetraquark states with variation of the Borel 
parameter M^. The A, B, C, and D denote the ccqq, cess, bbqq, and bbss channels 
respectively. The (I) and (II) denote the 6*7^ — 6*7'^ type and the 6*7^75 — 6*7^^75 
type interpolating currents respectively. 
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meson-meson type currents, there are contributions from the two-particle and many- 
particle reducible states, those contaminations are supposed to be small enough to 
be neglected safely. In fact, those contributions maybe considerable (and even out 
of control) and impair the predictive ability. In this article, we take the single pole 
approximation for the hadronic spectral densities, our predictions depend heavily on 
the two criteria (pole dominance and convergence of the operator product expansion) 
of the QCD sum rules; the numerical results are rather good. 
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Appendix 

The spectral densities at the level of the quark-gluon degrees of freedom: 
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where = ^ ^ ^ , = 2 ^ , A = m| = tuq = 

and A = 4(mQ + m,)^. 
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